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ABSTRACT 



This paper continues the analysis of a set of general relativistic 3D MHD 
simulations of accreting tori in the Kerr metric with different black hole spins. 
We focus on bound matter inside the initial pressure maximum, where the time- 
averaged motion of gas is inward and an accretion disk forms. We use the flows 
of mass, angular momentum, and energy in order to understand dynamics in this 
region. The sharp reduction in accretion rate with increasing black hole spin 
reported in Paper I of this series is explained by a strongly spin-dependent out- 
ward flux of angular momentum conveyed electromagnetically; when a/M > 0.9, 
this flux can be comparable to the inward angular momentum flux carried by the 
matter. In all cases, there is outward electromagnetic angular momentum flux 
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throughout the flow; in other words, contrary to the assertions of traditional ac- 
cretion disk theory, there is in general no "stress edge," no surface within which 
the stress is zero. The retardation of accretion in the inner disk by electromag- 
netic torques also alters the radial distribution of surface density, an effect that 
may have consequences for observable properties such as Compton reflection. 
The net accreted angular momentum is sufficiently depressed by electromagnetic 
effects that in the most rapidly-spinning black holes mass growth can lead to 
spindown. Spinning black holes also lose energy by Poynting ffux; this rate is 
also a strongly increasing function of black hole spin, rising to > 10% of the 
rest-mass accretion rate at very high spin. As the black hole spins faster, the 
path of the Poynting ffux changes from being predominantly within the accretion 
disk to predominantly within the funnel outffow. 

Subject headings: Black holes - magnetohydrodynamics - instabilities - stars: accretion 

1. Introduction 

In this paper, we continue the analysis of a scries of simulations of accretion disks in 
the Kerr metric introduced in Dc Villicrs, Hawlcy, & Krolik (2003; hereafter Paper 1). The 
emphasis in this paper will be on the properties of the accretion ffow inside the initial pressure 
maximum, where the net time-averaged motion is towards the black hole. This region of 
the accretion ffow, in the terminology introduced in Paper 1, comprises the inner part of 
the main disk body, the inner torus, the plunging region, and the corona. It is arguably 
the most important region of all in an accretion disk because it is here that most of the 
energy is released that ultimately powers the observed radiation and outffows. Our goal is to 
understand the peculiarities of dynamics in this region; our principal tool will be to analyze 
the ffows of three conserved quantities, mass, energy, and angular momentum. By so doing, 
we aim to shed light on the accretion process, how it is moderated by black hole spin, and 
how energy and angular momentum can be extracted from the spinning black hole. 

Accretion onto black holes differs from accretion onto ordinary objects in several ways, 
but two stand out: the absence of stable circular orbits inside of a critical radius, and the 
presence of frame-dragging due to black hole spin. The absence of stable circular orbits 
between the radius of marginal stability Vms and the event horizon means that in this region 
matter ffows in rapidly and the mass surface density is much smaller than in the main disk 
body where matter orbits stably. We may therefore talk about the disk having an "inner 
edge" in the vicinity of Vms- However, as emphasized by Krolik & Hawley (2002), where one 
places this inner edge depends on precisely which disk property is under consideration. For 
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example, two distinct edges can be defined by dynamical properties: the turbulence edge, 
where a transition takes place in the character of magnetic field evolution — from turbulence- 
driven evolution in the disk proper to evolution by magnetohydrodynamic (MHD) fiux- 
freezing, and the stress edge, within which the stress in the accreting fluid goes to zero so 
that there is no outward transport of angular momentum. Another two are more closely 
related to photon emission properties: the reflection edge, where production of Compton 
reflection features ceases, and the radiation edge, where intrinsic luminosity detected by 
distant observers goes to zero due to reduced radiative efficiency or relativistic effects. Fully 
relativistic simulations can explore the relation between these observable disk edges, and the 
detailed dynamics of the accretion flow. Defining the location and nature of these edges is 
central to linking disk dynamical models to disk observables because these edges determine, 
among other things, the global radiative efficiency of accretion. 

The other salient feature of relativistic accretion is the presence of frame-dragging when 
the black hole spins. Frame- dragging introduces a new dynamical effect into the accretion 
flow, and the spin of the black hole represents a potential source of energy output beyond 
that of the accretion flow itself. In a previous paper of this series (De Villicrs ct al. 2004) we 
presented evidence from the simulations that black hole spin energy helps to power axial jets. 
Here we focus on the flow of angular momentum and energy through the disk and corona, 
the influence of the black hole spin on the main accretion flow, and how the character of 
that flow depends on that spin. We flnd that frame-dragging and its associated effects can 
have a dramatic impact on the rate of accretion, the structure of the inner disk, and the 
fundamental flow of energy and angular momentum. 

This paper is organized as follows: After a quick summary of the simulations and the 
numerical diagnostics we used to study them (§2), we devote one section each to the flow 
through the disk of the three conserved quantities, mass (§3), angular momentum (§4), and 
energy (§5). In each case, we try to highlight how much goes where, how it is split between 
matter and electromagnetic contributions, and how the situation changes with changing black 
hole spin. Section 6 summarizes our flndings and comments on astrophysical implications of 
the numerical results. 



2. Preliminciries 

2.1. Simulation definition 

We solve the equations of ideal MHD in the metric of a rotating black hole. The specific 
form of the equations we solve, and the numerical algorithm incorporated into the GRMHD 
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code are described in detail in De Villiers & Hawley (2003). For reference we reiterate here 
the key terms and the definitions of the primary code variables. 

We work in the Kerr metric, expressed in Boyer-Lindquist coordinates, (i, r, ^, 0), for 

2 2 2 2 2 

which the line element has the form, ds — gu dt +2 gt^ dt d(f) + g^r dr + gee d9 + g^^ d(f) . 
We use the metric signature (—,+,+,+). The determinant of the 4-metric is g, and \/—g — 
q;^/7, where a is the lapse function, a = ^f**, and 7 is the determinant of the spatial 

3- metric. We follow the usual convention of using Greek characters to denote full space-time 
indices and Roman characters for purely spatial indices. We use geometrodynamic units 
where G = c = 1; time and distance are in units of the black hole mass, M. 

The state of the relativistic test fluid at each point in the spacetime is described by 
its density, p, specific internal energy, e, 4-velocity, If^, and isotropic pressure, P. The 
relativistic enthalpy ish = l + e + P/p. The pressure is related to p and e through the 
equation of state of an ideal gas, P = pe (F — 1), where F is the adiabatic exponent. For 
these simulations we take F = 5/3. The magnetic field of the fiuid is described by two 
sets of variables, the constrained transport magnetic field, Fj^ — [ijk] B\ and magnetic field 

4- vector V47r6'' = *F'^''U„. The ideal MHD condition requires C/''F^^ = 0. The magnetic 
field is included in the definition of the total four momentum, = {ph + ) WU^, 
where W is the Lorentz factor. We define the transport velocity as — [/*/[/*. 

In Paper I we presented results of a series of high- and low-resolution simulations, the 
KD (Keplerian Disk) set of disk models. These models have an initial condition consisting 
of an isolated gas torus orbiting near the black hole, with a pressure maximum at r ~ 25M, 
and a slightly sub-Keplerian initial distribution of angular momentum throughout. The 
initial magnetic field consists of loops of weak poloidal field lying along isodensity surfaces 
within the torus. In this paper the emphasis will be on the high-resolution models designated 
KDO, KDl, KDP, and KDE, which differ in the spin of the black hole around which they 
orbit, with a/M = 0, 0.5, 0.9 and 0.998 respectively These models used 192 x 192 x 64 
(r, 9, (j)) grid zones. The radial grid is set using a hyperbolic cosine function to maximize 
the resolution near the inner boundary, which is at Tin — 2.05 M, 1.90 M, 1.45 M, and 
1.175 M for models KDO, KDI, KDP, and KDE, respectively. The outer radial boundary 
is set to Tout = 120M in all cases. The ^-grid ranges over 0.045 tt < 9 < 0.955 tt, with 
an exponential grid spacing function that concentrates zones near the equator; reflecting 
boundary conditions are enforced in the ^-direction. The 0-grid spans the quarter plane, 
< < 7r/2. 
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2.2. Simulation diagnostics 

Three-dimensional numerical simulations generate an enormous amount of data, only a 
representative sample of which can be examined. Our analysis is based on a specific set of 
volume- and shell-averaged history data, taken every M in time, and complete data snapshots 
taken every 80M in time. Although the details of the history calculations are given in Paper 
I, we provide here a brief summary to clarify the calculations of mass, energy, and angular 
momentum transport in the inner torus and plunging region. 

The normalized shell-average of a quantity Q at radius r is calculated using 

{Q)A{r,t)^^j d9d(t>^gQ, (1) 

where the bounds of integration range over the full 9 and computational domains and A 
is the area of the shell. The density-weighted shell-average of a quantity will be denoted 
{Q)p = {p Q) / (p)- The flux JF of a given quantity through a shell at radius r is computed 
using a similar integral, but the result is not normalized to the shell area. For example, the 
accretion rate is defined 

M(r) = {pU') j de d(t)yf^pU' (2) 
In addition to the accretion rate (mass fiux) we consider the total energy fiux, 

{T\) = {ph U'- Ut) + {F\ Ft^) (3) 

4 TT 

and the total angular momentum fiux, 

{T\) = {phWU^)-r^{F\F^-) (4) 

4 TT 

The first term in each of these sums represents the contribution to the fiux from the fiuid. The 
second term represents the contribution due to the magnetic field, where the electromagnetic 
field strength tensor, F^^^ , reduces to combinations of transport velocities and constrained 
transport (CT) variables, as discussed in Paper I. 

We denote the time-average of a shell-averaged quantity (or flux) as 

{{Q))ir)-, dt{Q)ir,t). (5) 

''max ''min Jtmin 



The time-average of an integrated quantity (e.g., M) is denoted by a single bracket. The 
bounds of integration are chosen to analyze the late-time state of the accretion flow after the 
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initial transient described in Paper I has passed. Unless otherwise stated, tmin — 2000 M, 
which corresponds to about 2.5 orbits at the initial pressure maximum, and tmax = 8100 M, 
the end-time of the simulations. Such time-averages are important in establishing persistent 
features in various diagnostics; by their nature, these averages suppress short timescale 
variations, so we will also be interested in rms fluctuations about the time averages. 

The time-history data contains the total contribution on a shell, but docs not distinguish 
between contributions from bound material (i.e., in the disk and corona) and unbound (in 
the axial funnel and funnel- wall jet). For the purposes of this paper, where we are focused 
on the accretion flow, this is an important distinction. Thus, for many quantities, we will 
compute time-averages using the full data dumps, with integrations restricted to zones where 
the specific total energy satisfies —h Uf < 1, i.e., including only the contributions from bound 
material. 



3. Mass Accretion 

3.1. Inflow equilibrium and fluctuations about equilibrium 

The initial conditions in these simulations consist of an isolated torus of gas orbiting 
around the black hole. A relevant question is the degree to which the resulting accretion 
flow from such an initial state replicates a time-steady accretion disk. Outside the initial 
torus pressure maximum (r = 25 M), the gas must move outward as it accepts angular 
momentum from matter at smaller radius. This portion of the flow certainly cannot represent 
the behavior of a statistically time-steady accretion disk fed from large distances. Inside the 
initial pressure maximum, however, there is inflow throughout the simulation. Because the 
inner edge of the initial torus is at r = 15 M, well outside the marginally stable orbit, a 
significant amount of time is required before the flow reaches the black hole and a Keplerian 
disk forms. After that point in time we can investigate whether the accretion flow has 
reached an approximate statistical equilibrium. 

We first examine the rate at which gas is fed into the region r < 15M, namely the 
region inside of the inner boundary of the initial torus. Some sense of the time history can 
be seen in the spacetime diagrams of density (Fig. 5 of Paper I) which show both that density 
in the inner region is a highly dynamical quantity and that higher densities are found near 
the black hole in the high-spin models. Another view is given in Figure 1, which presents 
the total mass inside of r = 15 M as a function of time for each of the four models. In all 
cases significant matter infall begins by t = 500 M and the mass begins to build up. After 
t = 2000 M in the a/M = 0, 0.5, and 0.9 models, the disk in this region appears to settle 
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into a rough steady state of constant total mass. The high-spin model (a/M = 0.998) stands 
out: in that model, the mass inside r — 15 M continues to grow throughout most of the 
simulation. 

More detailed information can be obtained by examining the time-averaged accretion 
rate as a function of radius (M(r)). We begin the time average after t — 2000 M, i.e., after 
the first 25% of the total evolution time. Figure 2 plots (M(r)) for each of the four models. 
If the disk is in a statistically steady state, the time-averaged accretion rate should be nearly 
constant as a function of radius. This is clearly true inside of r = 15 M in the slowly-rotating 
simulations {a/M = 0, 0.5), and nearly so for the a/M = 0.9 model, but is not the case for 
the a/M = 0.998 simulation. Combining this result with Figure 1, we conclude that a state 
of local accretion equilibrium was attained for r < 10-15 M in three of the four simulations. 
Although the a/M — 0.998 simulation ran for the same amount of time, namely 10 orbits at 
r = 25M and hundreds of orbits in the inner disk, it did not reach a time-steady state with 
regard to mass flow. 

Despite the constancy with radius of the time-averaged accretion rate, the instantaneous 
accretion rate varies considerably with radius and with time. For example. Fig. 14 of Paper 
I shows that the accretion rate into the black hole varies substantially as a function of time, 
with significant power across a wide range of timescales. The dashed curves in Figure 2 show 
the rms deviation about the local mean accretion rate as a function of radius. The relative 
magnitude of accretion rate fluctuations (as compared to the mean) increases toward larger 
radii. Some of this effect may be attributed to the initial conditions and the absence of a 
mean accretion rate for r > 20M. However, a part of this effect is also likely to be real. It 
is significant, for example, that in each case the fluctuation level begins to rise well within 
the region where the accretion rate is nearly constant as a function of radius. 

There is also a trend for the relative size of accretion rate fluctuations at a flxed radius 
to increase with increasing black hole spin. Whereas Mrms/ (M) ~ 0.5 at r ~ lOM when 
a/M < 0.5, it rises to ~ 0.7 when a/M > 0.9; the trend is similar at other radii well inside 
the initial pressure maximum (r = 25M). Given that the mean accretion rate itself is poorly- 
defined when a/M — 0.998, it is not surprising that the rms fluctuation in the accretion rate 
in this simulation is comparable to the mean at all radii < 20M. 

Figure 2 also illustrates a point made in Paper I: the fraction of the initial torus mass 
accreted in a flxed time declines with increasing black hole spin. While 14% was accreted 
over the duration of the a/M — simulation, only 3.5% was accreted when a/M — 0.998. 
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3.2. 



Inflow time 



The mass accretion rate describes the bulk flow of matter toward the black hole; the 
inflow time is, on the other hand, a property of individual particles that make up that matter. 
As with many quantities in fluid dynamics, the inflow time may be studied from either the 
Lagrangian or Eulerian points of view. In this context, the "Lagrangian" inflow time is the 
(mean) time for fluid elements to move from radius r to the inner boundary; the "Eulerian" 
inflow time is the characteristic time for fluid elements to move past a flxed radius r. The 
former is of most interest for answering questions such as, "Does radiation have time to 
escape before a fluid cell passes the event horizon?" ; the latter is of most interest for gauging 
the time for the surface density at some radius to equilibrate with accretion fluctuations. 

We compute both inflow times from the mean inflow velocity obtained from a time- 
averaged density-weighted shell-integral. The quantity {{U^)p) is shown in Figure 3a for all 
four simulations. The inflow velocity is smaller with greater spin a/M at all radii less than 
about r = 15 M. The slope as a function of r also decreases with increasing black hole spin. 
Note too that, except for the highest-spin model, {{U'^)p) is very nearly a pure power- law in 
radius that shows no feature whatsoever at Vms- 

Because it is primarily used for comparisons to equilibration rates, it is most convenient 
to define the Eulerian infiow time in terms of the mean infiow rate. 



Figure 3b shows the radial dependence of this inflow time relative to the orbital period. In 
all cases but the most rapidly spinning black hole, this ratio rises gradually from ~ 2 near 
r = Vms to ~ 10-15 at large radii. The a/M — and 0.5 models are nearly identical; the 
a/M — 0.9 model shows a slightly flatter curve and a larger value inside Vms- The behavior 
of this ratio for the a/M — 0.998 model sharply contrasts with that of all the others: it is 
10-20 at all radn. 

For a Lagrangian interpretation, one can define the proper (i.e., fluid frame) radial infall 
time from a given radius to the inner radial boundary by integrating 



with ds^ — Qtt dt^ -\- Qrr dr^ using the time-averaged inflow velocity. Doing so produces curves 
that are very similar to those in Figure 3b, except that the inflow time is reduced by a factor 
~ 3-4. This contrast between the Eulerian and Lagrangian deflnitions is due to the sharp 
fall in the ratio of inflow time to orbital period at smaller radii. The inflow velocity scales 






(7) 
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roughly oc for the a/M — case, and is only slightly less steep when a/M — 0.5 or 0.9. 
Because of this rapid increase in velocity near the hole, most of the Lagrangian inflow time 
is spent at large radius. 

The Lagrangian inflow time provides one estimate for the location of the turbulence edge, 
the point where the magnetic field dynamics switch from being controlled by turbulence to 
being controlled by flux freezing in a plunging infall. Krolik & Hawley (2002) examined the 
ratio of infall to orbital time near the location of the turbulence edge, and found that it lies 
near the point where the (Lagrangian) infall time becomes shorter than the dynamical time. 
This occurs at ~ 1.3-1.4r„,,,5 for all three models with a/M < 0.9. When the black hole spins 
very rapidly, however, this ratio is almost constant as a function of radius, not betraying a 
clear location for the turbulence edge at all. 

Indeed, a striking aspect of Figure 3b is the contrast between the three more slowly- 
spinning cases and the most rapidly-spinning one. The slow-spin cases behave as one would 
expect: well outside the marginally stable orbit, inflow takes multiple orbits because it is 
limited by the slow loss of angular momentum due to magnetic torques. As one approaches 
the location of the marginally stable orbit, however, there is less to prevent inward radial 
motion at constant angular momentum because the effective potential is very nearly flat. 
There, the ratio of infall time to orbital period drops toward one. 

The a/M — 0.998 case is very different. This contrast presents a puzzle because the 
effective potential seen by its material is not qualitatively different from the others. As was 
discussed in Paper I (see Fig. 9) , the specific angular momentum distribution in the accretion 
disk is close to Keplerian in all four simulations. It is therefore difficult on the basis of 
gravitational orbital dynamics alone to understand why the ratio of inflow to orbital time 
remains large even at r = r^^ when a/M = 0.998, and why the transition from turbulence 
to plunging inflow is inhibited. We believe that this result has implications for the impact 
of electromagnetic stresses on the angular momentum budget. This will be discussed in § 4. 

3.3. Radial distribution of surface density 

One deflnition of an accretion disk's inner edge is a rapid decline in the disk's surface 
density distribution. Figure 4 is the time-averaged rest-mass surface density distribution 
as a function of radius for all four simulations. The mean surface density for each of the 
simulations is normalized to the initial torus mass for that simulation. In addition, in order 
to avoid creating a misleading visual impression, the radius in this plot is proper radius, 
not Boyer-Lindquist coordinate radius; the relation between this radial variable and Boyer- 
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Lindquist radius is defined as 

= / dr' V9rr{r'), (8) 

Jrir, 

where the metric coefficient is evaluated in the equatorial plane. We derive the surface den- 
sity from the history data, which provides the total mass in each radial shell for every M 
in time. To obtain a surface density we divide the shell mass by an equatorial unit area 
and average over time for t > 2000 M. We restrict radial coverage to regions within 15M 
in order to consider only those places in approximate inflow equilibrium; for this reason, 
when plotted relative to Rms, the curves extend for different distances. All these curves 
may be compared with one illustrating a conventional disk model (triple dot-dashed line). 
This curve was constructed combining the run of stress predicted by the (zero-stress bound- 
ary condition) Novikov-Thorne model and the fixed ratio of integrated stress to integrated 
pressure suggested by Shakura & Sunyaev (1973). 

Several points stand out in Figure 4. All four simulations show quite similar behavior, 
with S rising steadily with increasing radius, steeply within the plunging region and more 
gradually in the disk proper. As was emphasized by Krolik & Hawley (2002), there is no 
sharp break in surface density, not even at the marginally stable radius. Normahzed to the 
initial mass, the surface density in the marginally stable region declines by about a factor 
of two as a/M increases from to 0.9, but then rises again when a/M — 0.998. In sharp 
contrast to this behavior, the conventional model predicts a well-defined maximum in E 
at 2-3 Rms, with a drastic drop-off at and inside and a more gradual decline at larger 
radius. The reason for the very sharp break in the standard model is that the net radial 
accretion velocity is much smaller than either the orbital velocity or sound speed outside 
Rms, but is assumed to abruptly rise to the free-fall velocity inside Rms- The result is an 
abrupt drop in surface density. As can be seen in Figure 3a, the accretion velocity in our 
simulations increases smoothly and continuously from outside Rms inward toward the hole. 

The surface density distribution is important in locating the reflection edge of the disk, 
the point within which the disk can no longer efficiently Compton scatter X-rays (Krolik 
& Hawley 2002). Compton reflection is a process whereby dense, cool matter in the disk 
reflects hard X-rays generated somewhere outside the disk, imprinting upon their spectrum 
the signature of absorption by heavy element K-edges. Efficient reflection requires a disk 
column density large enough to make it thick to Compton scattering. Because the column 
density must decrease sharply with diminishing radius in the plunging region, the reflection 
edge is generally assumed to be located at or near r^s- 

Although the simulations have no absolute density scale, the results may be phrased 
in terms of an absolute density once the accretion rate in Eddington units is specified: the 
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relation is E = 38 m g cm~^, where the accretion rate relative to Eddington m is defined 
assuming unit radiative efficiency in rest-mass units {MecLcl = 1-7 x 10^^ mM/M© g s~^). 
A maximum value for fn in these units would be ~ 20 (if the efficiency is as low as 0.05 
and the luminosity is exactly Eddington), and in most instances one would expect it be 
a good deal less. The surface density in code units for the three slower-spinning simula- 
tions near Vms is — 0.01, while it is ~ 0.06 for the a/M = 0.998 simulation. Translating 
to physical units, we predict Compton optical depths tt — 1.3(m/10)(S/0.01 code units). 
Thus, rapidly-accreting systems may indeed find their reflection edges near r^s, but in more 
slowly-accreting systems it may be significantly farther out. And, of course, in all cases, 
the location of this edge can vary significantly in time. Because disk thermodynamics can 
infiuence the radial distribution of surface density, it may be too early to say definitively 
how the efi^ects we are examining alter classical expectations, but it is clear that they can 
strongly affect this basic property of disks. 

To summarize the results of this section, we find that the low spin models with a/M < 
0.9 develop accretion disks that are in a quasi-steady state, but the highest spin model, 
a/M = 0.998 does not. For fixed surface density at large radius, the accretion rate into the 
hole decreases with increasing a/M, although the accretion rates feeding mass into the region 
inside the initial torus are comparable. This means that the mass of the inner torus tends 
to be greater for larger black hole spin. Accretion infall times, too, are longer with larger 
spin. All these results suggest that the spin of the black hole affects the disk in ways that go 
beyond simply determining the location of the marginally stable orbit. The mechanism by 
which the black hole does so is elucidated by the next topic, the angular momentum fluxes 
in the disk. 



4. Angular Momentum Transport 

The stress tensor 

Ti: = phU^U, - + WhW^'U^'U, + (p+ ||6||V2)< (9) 

describes the flux of four-momentum. Conservation of four-momentum means that V ^Tj^ = 
0, where the V indicates covariant differentiation. The four terms contributing to the stress 
tensor each have simple interpretations: the first is the flux directly associated with matter; 
the second is the purely electromagnetic part due to correlations in the magnetic fleld; the 
third is the part proportional to the magnetic energy advected with the flow; the fourth is 
the momentum flux associated with pressure. In the context of accretion disks, the most 
interesting component of this tensor is T^, the angular momentum flux in the radial direc- 
tion. The second term then represents the magnetic torque that moves angular momentum 
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outward, while the first and third represent the angular momentum moving inward with the 
accretion flow, and the fourth is zero. Thus, the magnetic torque contribution —}fh^ has 
particular importance for accretion disks because it is responsible for transferring angular 
momentum outward from one fluid element to another. The net flux (the sum of all the 
terms) is the conserved angular momentum flux through the disk. 



To understand better the flow of angular momentum in the disk, it is helpful to look at 
the spherical shell-integrated time-average of each of the three terms in equation 9 separately 
as a function of radius (Fig. 5). Each of the three terms is integrated over the angular 
coordinates, but we restrict the integration to the region where the flow is bound, i.e., 
—hUt < 1. This removes the contributions from the magnetic fields in the funnel which, as 
discussed in De ViUiers et al. (2004), are not insignificant when the black hole is rotating. 
The computation is done using complete data sets spaced every 80M in time; these are then 
averaged over the last 75% of the simulation. 

Because V^T^^ = 0, the integral over a radial shell, J dflT^, is constant in a steady-state 
disk. Separating the terms of the r-0 component of the stress tensor and employing Gauss's 
Theorem transforms this statement into 



In conventional disk models it is assumed that the stress term in equation (10) is zero inside 
some radius that lies well outside the event horizon. Generally this "stress edge" is assumed 
to lie at the marginally stable orbit, or, in more detailed hydrodynamical models, it is taken 
to be the sonic point, which is generally a short way inside r^g (e.g., Abramowicz et al. 1988). 
However, as Page & Thorne (1974) remarked, and Krolik & Hawley (2002) discussed in the 
context of pseudo-Newtonian models, when magnetic fields are important, the stress edge 
need not coincide with rms- We will shortly examine where it may be found in our general 
relativistic simulations. 

Time-steadiness also implies that the accretion rate is constant as a function of radius, 
so equation (10) requires the shell-integrated stress to increase outward to match the angular 
momentum per unit mass, which is oc r^^^ in the Newtonian limit. Physically, this means that 
for matter to move inward it must lose angular momentum by transferring the appropriate 
amount of angular momentum outward. A steady state at a fixed position is maintained 
when that transfer is balanced by matter fiowing in carrying angular momentum. 



4.1. Shell-integrated angular momentum flux 




(10) 
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We can directly compare this picture to the simulation results. In Figure 5 we see that 
in the low spin cases, a/M = 0.0 and a/M = 0.5, some of this picture is reproduced in a 
time-averaged sense. The shell-integrated stress rises outside the marginally stable orbit, 
accompanied by a compensating gradient in the matter flux. However, rather than going to 
zero inside r^s, the shell- integrated stress, while smaller than in the main disk body, stays 
at a constant level. 

In the two higher spin models the departures from conventional expectations are even 
greater. Just as in the low-spin cases, stress continues throughout the plunging region, but 
its magnitude (there and in the disk proper) is everywhere proportionally larger. In rough 
terms, its magnitude is comparable to the angular momentum flux carried with the matter. 
We conclude, therefore, that, as viewed in the coordinate frame, there is no stress edge; 
stress continues all the way from the disk body to the event horizon. 

A second, even more dramatic, contrast with conventional views is that in the inner 
regions the gradient in the shell-integrated stress has the wrong sign: the magnetic angular 
momentum flow decreases outward. Where this is the case, rotationally-supported matter 
does not on average lose angular momentum by magnetic torques. This change in sign of the 
gradient of the angular momentum flow appears to be the reason why the rate of accretion 
onto the black hole is smaller in the two simulations with rapidly-spinning black holes, and 
likewise why the mass of the inner torus becomes large in these two cases. 

4.2. Mean specific angular momentum 

Additional information can be obtained through an analysis of the mean accreted angu- 
lar momentum per unit rest-mass, j. Again with the integrals restricted to the bound matter, 
this is deflned as j(r) = ((T^))/(M). We plot the function j(r) in Figure 6. Here the solid 
line shows the specific angular momentum associated with the total angular momentum fiux, 
while the dashed line is that of the matter fiux alone. The a/M = and a/M = 0.5 simu- 
lations behave very much in line with expectations; j is very close to constant with radius 
through the disk, indicating that the angular momentum fiux has equilibrated to an approx- 
imately time-steady state throughout these regions. Although some angular momentum is 
carried in the jet when a/M — 0.5, it is too little to have much effect on the total. In each 
model the dotted line indicates the value j{rms) associated with the marginally stable orbit, 
and it is clear that jin < j{rms), where jin is the value carried into the black hole. This 
demonstrates that torques inside the marginally stable orbit do remove some angular mo- 
mentum from material fiowing through that region, but the effect is not dramatically large: 
it ranges from 6% to 8%. 
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In the two simulations with the more rapidly-spinning black holes, the radial slope in 
j{r) provides evidence that they are not in steady state, particularly the a/M = 0.998 
model. Again the specific angular momentum carried by the mass does not differ much 
from the value j(r). However, in these cases there is a significant difference between the net 
angular momentum carried by the bound material and the total. The angular momentum 
flux associated with the jet is ~ 10% of that carried in bound matter when a/M = 0.9, but 
is comparable to the flux through the bound matter when a/M — 0.998. When all is said 
and done, j is 12% below j{rms) for the a/M = 0.9 model, and 42% below for a/M = 0.998. 

If the conditions of accretion remain constant long enough for the accreted mass to 
dominate the original mass in the black hole, j determines a/M. As noted in Paper 1, with 
jin = 0.8 in the a/M — 0.998 simulation, in the long run the spin of the black hole can not 
be maintained (see also Gammie, Shapiro, & McKinney 2004). 

4.3. Angular distribution of angular momentum flux 

We can see more precisely where the angular momentum goes by studying the angular 
distribution of the magnetic torque contribution to the stress, T^{r,9). In Figure 7 the 
dashed lines are the (inbound) matter component, plotted at 2r^s and Vms, and the solid 
lines are the (outbound) total electromagnetic component, plotted at the same locations as 
the matter component and at the inner boundary as well. 

The plot for the Schwarzschild model is straightforward. The matter angular momen- 
tum flux narrows going from 2r^s to r^^, while the magnetic stress peaks through the disk 
and drops off through the corona. In the corona there is a region where the outward elec- 
tromagnetic stress exceeds the inbound matter flux. There is neghgible electromagnetic 
angular momentum flux at the inner boundary, and what there is is ingoing, represented by 
a dot-dashed fine. 

In contrast, all the non-zero spin simulations have a signiflcant outward angular mo- 
mentum flux carried by magnetic flelds throughout the region between the marginally stable 
orbit and the inner boundary. Moreover, as was also shown in Figure 5, the magnitude of 
this flux relative to the stress associated directly with matter {phlf^U^, the dashed curve 
in the figure) increases steadily with increasing a/M. As a/M increases, the angular width 
through which the magnetic angular momentum fiux travels both spreads and develops a 
double-maximum structure. That is, when the black hole spins slowly, the magnetic angular 
momentum flux is largely conflned to the region where the matter flux is greatest, a span of 
perhaps ±30° on either side of the equatorial plane. Moving from the disk to the corona. 
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there is a point where the magnitude of the magnetic stress exceeds that of the matter stress. 
When the black hole spins rapidly, the magnetic angular momentum flux is relatively con- 
stant with polar angle through the disk region, but rises to a maximum at the funnel wall 
and slowly declines in the funnel interior toward the axis. 

In the Schwarzschild simulation the magnetic stress peaks within the disk, and falls as 
the flow approaches the hole, becoming inward- directed at the inner boundary. In the a/M = 
0.5 simulation the magnetic stress can be described in terms of two pieces: a contribution 
spanning a broad angular range that increases slowly with increasing radius, and another 
contribution whose angular extent roughly matches that of the disk matter and grows toward 
larger radii. As the black hole spins faster, the relative importance of the broad angular range 
component grows, and a double-peak structure develops, with the greatest stress occurring 
near the funnel wall. 

We can understand these trends in terms of a changing balance between effects arising 
from the magnetorotational instability in the disk and effects stemming from magnetic cou- 
phng to the spin of the black hole itself. In the disk, the stress results from MHD turbulence 
generated by the MRI (Balbus & Hawley 1998), and the magnetic stress is tied directly to 
the presence of the matter. The narrow disk component can be identified with this process; 
it increases with radius because its gradient must match that of Mf/^, which increases with 
radius. The broader component we associate with frame- dragging acting on the magnetic 
field in the corona and outflow. Close to the black hole, the magnetic fleld intensity varies 
little with angle, so there can be stress across a wide range of polar angles. Its magnitude 
changes little with radius because the matter density is so small in the corona and outflow 
that there can be little angular momentum exchange between fields and matter. 

In rough terms, it is possible to separate the electromagnetic angular momentum flux 
going into the disk from that going into the corona by evaluating this flux on the shell at 
^ = ^ms and distinguishing disk from corona in terms of polar angle. Applying this criterion 
quantifies some of the remarks of the previous paragraphs: as a/M increases, the ratio of the 
electromagnetic angular momentum fiux entering the corona to that entering the disk rises 
from tiny (< 10~^ in the Schwarzschild case) to ~ 0.3 (a/M = 0.5), ~ 1 (a/M = 0.9) and 
~ 3 {a/M = 0.998). Because angular momentum can readily move from matter to fields and 
back again, and because polar angle distribution at a fixed radius does not fully distinguish 
the ultimate end of the fiux, these numbers are not entirely well-defined; nonetheless, they 
do express the qualitative point that the corona receives a growing proportion as the black 
hole spins faster. 

When the black hole spins rapidly, the overall level of electromagnetic angular momen- 
tum fiux is so high that, despite the growing share devoted to the corona, the portion going 
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into the disk is large enough that it fundamentally alters the angular momentum budget of 
the accretion flow. It appears that in both the a/M = 0.9 and a/M — 0.998 simulations, the 
time-average magnetic torque adds to the matter's angular momentum over a substantial 
radial range from near the hole to r^s and beyond (Fig. 5). As a result of this process, the 
accretion flow is retarded. In Paper 1, we observed that in each simulation an inner torus is 
created in the region immediately outside Tmsi and that the prominence of this inner torus 
increases with increasing a/M. The electromagnetic stress driven by the spin of the black 
hole is the likely origin of this dependence on a/M. In the case of the a/M — 0.998 sim- 
ulation, magnetic torques are so strong that the mean of the local angular momentum can 
exceed the value associated with a circular orbit. For example, the density-weighted ratio of 
local angular momentum to the angular momentum of a circular orbit at that radius rises 
from ~ 0.95 to ~ 1.05 in the inner torus region (see also Fig. 13 of Paper I). The effective 
potential for the flow is then not as flat as expected in the marginally stable region, and 
inflow cannot occur as easily (see also Fig. 3). 



4.4. Fluid frame stresses 

It is also useful to look at the stresses in the fluid frame. The MHD turbulence is 
primarily due to local instabilities, so its physics is best analyzed in the local fluid frame. In 
addition, Novikov & Thorne (1973) showed that, if the disk is assumed to be time-steady, 
axisymmetric, and geometrically thin, the vertically-integrated stress in the fluid frame can 
be written in the physically-enlightening form 

W;{r)^^RT{r/M,a/Ml (11) 

where the function Rt encapsulates both relativistic corrections and the net angular momen- 
tum flux (the notation here follows Krohk 1999). The assumption that stress goes to zero 
at the marginally stable orbit means that Rt{i^) = at all radii r < r^s- In the Newtonian 
limit of large r, R^ — > 1. In this relation, "stress" means the part of the stress tensor other 
than that directly associated with accretion, i.e., — — MUii,{r), where C/^(r) is the 
circular orbit angular momentum at r. 

To compare our results with the Novikov- Thorne prediction, we compute a closely- 
related quantity. 



r ded(f)y/^[\\b\\^U^'U^-lfK] e^e'f. 

^I2ir) = r \ (12) 
! d4>^-g{e^7:/2)/gee{e^7:/2) 



where the e^^^ are an orthonormal tetrad designed to represent the directions labeled by {v) in 
the fluid frame (see Appendix for details about their construction) . We limit the integration 
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to the bound portion of the flow, where —hUt < 1. We divide the shell-integrated fluid-frame 
stress by the area per unit radial coordinate in the equatorial plane in order to give this 
the same units as W^. 

We compute W^^j (r) at t — 8000 M for all four simulations; this is shown in Figure 8. 

To facihtate comparison of simulations with different black hole spins, the value of W^^j dis- 
played in each panel is normalized to the initial maximum pressure found in that simulation. 
There is rough agreement with the Novikov-Thorne prediction (assuming the time-average 
accretion rate at the inner boundary) for LSr^s i$ ^ i$ 40M. At larger radii, one would not 
expect agreement because that is the region in which the matter in our simulations must 
move outward to absorb the angular momentum flux, whereas in the Novikov-Thorne model 
that region is supplied by accreting matter from farther away. Inside r^s, however, where 
the Novikov-Thorne model predicts zero stress, our data are dramatically different. We find 
that the stress continues quite smoothly, and, in fact, continues to increase in amplitude 
from to the inner boundary whenever a/M > 0. Only in the Schwarzschild case does 
the fluid-frame stress eventually go to zero, albeit only very close to the event horizon. 

In other words, there is considerable outward angular momentum flux in the fluid frame 
virtually throughout the flow, including the plunging region. Moreover, whenever a/M > 0, 
there is always at least some outward angular momentum flux in the fluid frame even at the 
smallest radii outside the event horizon. Further, even with the overall normalization, the 
general level of the fluid- frame magnetic stress rises substantially with increasing a/M. This 
is due to the increasing amplitude of the magnetic fields. Thus, while there is a fluid-frame 
stress edge in the Schwarzschild case just outside the event horizon, there is no fluid-frame 
stress edge when the black hole spins. 

The standard Novikov-Thorne model assumes a cold, unmagnetized disk, with h = 
= 0. The origin of the stress is unspecified, although it is presumed due to turbulent 
motions. In our simulations neither h nor is zero. In the calculation of Figure 8 we 

included only the magnetic stress and neglected the turbulent Reynolds stress, i.e., phU^5U^, 
where the perturbation is with respect to a local circular orbit. Local simulations (Hawley, 
Balbus & Winters 1999) show that the Reynolds stress is generally a factor of several less 
than those due to the magnetic fields. Obviously the presence of magnetic fields explains the 
quahtatively different behavior of W^^j and from the marginally stable region inward. 

Another potential infiuence on the stress edge is the temperature in the disk and the 
associated vertical scale height. Purely hydro dynamical models of accretion disks more de- 
tailed than the Novikov-Thorne model (e.g., Muchotrzeb & Paczyhski 1982; Matsumoto et 
al. 1984) relax the assumption olh — Q and consider the effects of finite geometric disk thick- 
ness. These studies indicate that, although the inner boundary for stress may not coincide 
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exactly with the marginally stable orbit, in thin disks it should he very close. Even in "shm 
disks" (Abramowicz ct al. 1988), the stress does not extend all the way across the plunging 
region: in the pseudo-Newtonian approximation, designed to mock up Schwarzschild metric 
dynamics, the innermost radius for the end of stress stretches inward from 6M to 4M as the 
disk thickens. But it is clear that the purely hydrodynamic picture is inadequate; magnetic 
stresses operate in ways that Reynolds stresses and simple viscosity cannot. Further, in 
the MHD simulations there is significant magnetic field in the surrounding disk corona; the 
magnetic scale height is greater than the gas pressure scale height of the disk. In fact, the 
magnetic scale height near the black hole is ~ r more or less independent of the matter scale 
height in the disk. Thus, when magnetic fields are included, the picture changes qualitatively 
from the purely hydrodynamic predictions. 

Earlier we found that, when stress is measured in the coordinate frame, there is no stress 
edge anywhere. When it is evaluated in the fiuid frame, the situation changes only slightly: 
when a/M > 0, there is similarly no stress edge; when a/M — 0, one exists, but it is deep 
in the plunging region, immediately outside the event horizon. 

5. Energy transport 

5.1. Energy fluxes 

The conventional approach to defining the radiative efficiency of accretion onto black 
holes was set out by Novikov & Thorne (1973). In the standard picture, thermal energy 
in the disk is assumed negligible compared to orbital energy, in part because the thermal 
loss time is assumed to be very short compared to the inffow time everywhere outside Vms- 
With those assumptions, the radiation rate in each ring is simply the difference between the 
energy brought into that ring by mechanical work (torques acting on the ring by neighboring 
zones) and energy lost from that ring by accretion. Both contributions amount, of course, 
to the radial divergence of an energy ffux which does not vanish; the energy deposited by 
the ffux is the disk's luminosity. As the ffuid drops toward the hole, it becomes increasingly 
bound, with a binding energy equal to that of the circular orbit at its radius. The energy 
of the last circular orbit at r„is sets the canonical efficiency of standard thin disk accretion, 
rj = 1 + Utij-rns)- For a Schwarzschild hole this is 0.057. 

In the simulations, things are less straightforward. As in the standard picture, energy 
is extracted from orbital motion as the gas accretes, but that energy is not radiated. It first 
goes into magnetic fields and velocity fiuctuations within the turbulent disk. Some magnetic 
field rises out of the disk into the corona, and a general coronal outfiow removes some energy 
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from the disk. As the turbulence in the disk is dissipated at the grid scale, some of the 
energy is converted into heat by artificial viscosity; some is lost. Enough of this energy is 
retained to keep the disk hot. 

We can investigate how energy is transported in these simulations by considering the 
magnitude of the energy fiux, given by the component of the stress tensor. Just as for 
the angular momentum fiux, it has both a fiuid part 

(^/■)fluid = PhU'^Ut (13) 

and an electromagnetic part 

TOem = \\b\\'U^Ut-b^bt 

= g"' [f^S'l'B^ - g^^8^B'l' - g^'^'S'^S^) , (14) 

where £^ = V'^B'^ - V^B''', E'^ = V'B'^ - V'^B'', and £^ = V^B'' - V'B^. The last form for 
the electromagnetic energy flux shows that it is the radial component of the Poynting flux, 
modifled by a frame-dragging correction. 

In a steady-state disk, T/' integrated over spherical shells should be constant with radius. 
However, because its largest contribution is from the flow of rest-mass, its radial behavior 
is dominated by that flux, a conserved quantity to which we have already devoted much 
attention (§3.1). In order to focus attention on the non-rcst-mass portion, we instead study 
shell-integrals of time-averaged - M = {{p{hUt + l)t/^ + WbW^WUt - ^h)). In a time- 
steady disk, the constancy of M with radius means that this "adjusted" energy fiux should 
likewise be constant as a function of r. In an effort to make still closer contact with ordinary 
intuition, in Figure 9 we plot 1 — {{T[))/ (M), the binding energy per unit rest-mass accreted, 
but restricting the integral to cells in which the matter is bound. In this figure, the dashed 
line is the speciflc energy carried by the fluid term alone (eq. 13), and the solid line is 
obtained from both the fluid and electromagnetic components. As references, there are two 
dotted lines: the straight one shows the binding energy of the marginally stable orbit, the 
curved one the binding energy of a circular orbit in the equatorial plane at that radius. 

If all the simulations achieved a steady state within the region shown, and if there were 
no energy exchange between bound matter and anything else, the solid line would be flat. 
As is immediately apparent, it is not. At large radii, the fluid binding energy flux (dashed 
line) closely tracks the orbital energy curve. This is because the numerical energy loss time 
is shorter than the inflow time, crudely mimicking the effects of radiation losses in a real 
accretion disk. At smaller radii, the fiuid binding energy is smaller than the orbital binding 
energy for two reasons: artificial viscosity adds heat to the gas, and the infiow time becomes 



-20- 



shorter than the numerical loss time. Some energy is lost to the outflow, especially when 
a/M > 0.9, but not enough to strongly cool the disk. Here the mimicry is closer in spirit to 
a non-radiative disk. 



5.2. Poynting flux 

At all radii there is a substantial offset between the total energy flux curve and the 
fluid-only curve. This offset is due to the outward energy flux conveyed electromagnetically 
(Fig. 10). With the exception of the inner boundary in the zero-spin simulation (where it is 
also very small in magnitude), the Poynting flux is always directed outward. Independent 
of the black hole spin, there is always Poynting flux associated with the work done by 
torques in the disk. In the outer portions of the accreting disk (r ~ 15M), where our initial 
conditions make all four simulations very similar, its magnitude varies little with black hole 
spin. However, in the inner disk, where the mean magnetic fleld strength grows substantially 
with more rapid rotation, so does the Poynting flux. At the equator on the inner boundary, 
it increases by about a factor of two from a/M = 0.5toa/M = 0.9 and by almost an order 
of magnitude from a/M = 0.9 to a/M = 0.998. 

Poynting flux in the jet outflow varies far more dramatically with black hole spin. It is 
neghgible in the Schwarzschild simulation, while in the a/M — 0.998 simulation the Poynting 
flux in the funnel wall completely dominates the Poynting flux in the disk. Another way to 
make this comparison is to consider the ratio of the total (outward) shell-integrated Poynting 
flux to the (inward) fluid energy flux (nearly all rest-mass energy). This ratio is plotted in 
Figure 11. Here we see that relative to the rate at which fluid energy is accreted, the Poynting 
flux integrated over the inner boundary rises from 0.25% when a/M = 0.5 to ~ 1% when 
a/M = 0.9 to ~ 10% when a/M = 0.998. These ratios increase steadily outward for the 
a/M — and 0.5 simulations, and have a local maximum around r = 3-5 M for the high spin 
models. In the simulation with the most rapid spin at the radius where this ratio reaches its 
maximum, Poynting flux carries outward ~ 60% of the ingoing fluid energy flux. 

Another way to look at these results is from the point of view of the Penrose process. 
In its simplest form, this process depends on the fact that inside ergospheres there can 
exist negative energy orbits, so that black holes can give up energy (i.e., lose mass) by 
absorbing particles with negative energy. An analogous process exists for magnetic flelds. 
The electromagnetic "energy density at inflnity" may be deflned (Koide 2003) by 

e^u/a = -Tl{EM) = - (1/2)^" [g^^{rf + ^{8^ + g^^{£'^f] + (1/2) {g'^S^Y 

+ [g^^g^^iB'-f + /V^(i3^)2 + g'^-g^^B^f] . (15) 



(16) 



/,From this form we see that, if all the spatial elements of the inverse metric were positive 
definite, could never be negative. However, inside the ergosphere, g'^'^ < 0. Thus, 
negative electromagnetic energy-at-infinity is associated with regions inside the ergosphere 
where the components of the field tensor associated with poloidal magnetic field {B'^ and B^) 
and toroidal electric field {£'!') are large in magnitude compared to the other field components. 
Moreover, because S'^ — V^B^ — V^B^ , a large poloidal magnetic field also helps to make £^ 
large, particularly when the poloidal transport speed is substantial. That an electromagnetic 
version of the Penrose process is at work is suggested by the fact that in the a/M = 0.9 and 
0.998 simulations, the time-averaged e^-^ < in the region within 40° - 60° of the equatorial 
plane and ^ 0.3M of the inner boundary. 

One of the key questions of black hole accretion is whether or not energy can be extracted 
from the region inside r^g. Although it is difficult to say how much radiation will escape 
from inflowing gas within the plunging region (discussed below), it is clear that the Poynting 
flux provides an avenue for the outward transport of energy from the plunging region and 
from the black hole itself. 



5.3. Radiation 

Although the present simulations do not include radiative losses, it is interesting never- 
theless to consider whether the standard assumptions of thin disk theory would apply. The 
issues involve where and how rapidly energy is dissipated into heat, and whether there is 
sufficient time for the heat to be transformed into photons and radiated prior to the gas 
entering the black hole. Well outside the marginally stable orbit, the time required for these 
processes to occur should be significantly shorter than the infiow time. Where that is true, 
the conservation law approach to computing the radiation rate should work well. In the sim- 
ulations the turbulent energy released by dissipation is captured only partially by artificial 
viscosity; numerical losses thus mimic radiation losses, and in most of the main disk body, 
the numerical loss time is appropriately shorter than The only difference between the 
simulations and the Novikov-Thornc model in this regard is that we do not force the stress 
to go to zero at Vms- A formalism for tracing the consequences of that changed boundary 
condition in time-steady disks was developed in Agol & Krolik (2000). Applying it to the 
simulation results indicates an enhanced radiative efficiency from this part of the disk that 
is relatively modest for the cases with a/M — 0, 0.5, and 0.9. Because the a/M — 0.998 
simulation did not achieve an inflow steady-state, the assumptions of this formalism may 
not be valid, but the data suggest that the enhancement to the radiative efficiency may be 
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considerably larger. 

Inside ~ 2r^s, the inflow time falls rapidly (see Fig. 3), becoming roughly an orbital 
period near and inside r^g- K then becomes questionable whether the several processes that 
stand between inter-ring work and radiation have enough time to run to completion. If not, 
the gas advects its energy inward, much as envisioned by Rees et al. (1982), Narayan & Yi 
(1995), and Abramowicz et al. (1988), albeit for different reasons. However, particularly 
in the plunging region proper, new processes can lead to radiation. Machida & Matsumoto 

(2003) suggest that field reconnection in the plunging region could lead to X-ray fiares. 
Rather than magnetic field being dissipated on very small lengthscalcs at the end of a 
turbulent cascade, it can instead by dissipated in large-scale current sheets. In Hirose et al. 

(2004) we drew attention to the prevalence of high current density regions in the inner disk 
and plunging region where the current density is so high that rapid field dissipation might 
be expected. If the local temperature in these current sheets can be raised high enough that 
Compton scattering is the dominant cooling mechanism, radiation can release energy very 
rapidly: the Compton cooling time is 



where jj,^ ~ nip is the mass per electron and L is the luminosity of seed photons injected 
within radius r. Some of these photons will be captured by the black hole, and some will 
strike the disk farther out, but some portion will also reach distant observers. To compute 
the actual radiative efficiency of flows of this sort will require calculations that explicitly 
balance the various competing processes and timescales. 



In Paper I of this series we reported that the fraction of the initial mass accreted in 
a flxed time diminishes sharply with increasing black hole spin. We also noted that the 
angular momentum per unit rest-mass that is accreted likewise drops as the black hole spins 
more rapidly, and by signiflcantly more than the associated decrease in the speciflc angular 
momentum of the innermost stable circular orbit. On the basis of the more detailed analysis 
presented here, we can now explain both these effects. Not surprisingly, they are closely 
related. 

An examination of the angular momentum fluxes in the accretion flow reveals how 
significant magnetic fields are to the evolution of the inner region of a black hole accretion 
disk. The low a/M models attain a quasi-steady state that behaves much as one would 
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6. Discussion and Summary 
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expect from simple accretion disk theory. The net flux of angular momentum through the 
disk is essentially equal to the value carried into the black hole. Within the flow itself 
outward-going magnetic angular momentum flux balances the excess in the inward-going 
component carried by the matter flow. 

Contrary to standard accretion theory, however, there are always magnetic stresses 
exerted across the location of the marginally stable orbit and in the plunging region interior 
to this point. In fact, when the black hole rotates, the r-0 component of the stress tensor 
in the fluid frame rises more steeply with decreasing radius in the plunging region than in 
the disk itself. Thus, we have shown that the concerns voiced quietly in Page & Thorne 
(1974) and Thorne (1974) have been borne out: the presence of magnetic fields invalidates 
the "zero-stress" inner boundary condition. 

Hirose et al. (2004) showed that the relative magnitude of the magnetic field in the inner 
disk grows steadily with increasing black hole spin. This is a by-product of the creation of 
the dense inner torus: higher pressure gas can support stronger magnetic field when the 
field's origin lies in the magneto-rotational instability. Thus, as the spin increases, the 
magnetic stresses grow rapidly stronger. There is more going on than can be attributed 
merely to stronger field strengths, however. The rotating hole itself represents a source of 
angular momentum and energy that can be carried outward by the magnetic stress. When 
a/M > 0.9, the magnetic stresses carry out of the ergosphere an angular momentum flux that 
is comparable to what is brought inward with accreting matter. A striking suppression of 
the accretion rate is the result. In the most extreme case (a/M = 0.998), accretion through 
the inner disk is retarded to such a degree that even in a simulation time of 8100 M no 
accretion equilibrium is achieved; rather, the mass of the inner torus builds steadily through 
the simulation as matter accretes from outside, but is held back not far outside r^g- This 
result is reminiscent of the suggestion by van Putten and Ostriker (2001) that magnetic 
torques driven by rotating black holes might lead to a state of suspended accretion, but 
the fleld structure we see (turbulent field, largely confined to the accreting matter) is quite 
different from the one posited by them (organized poloidal field with loops spanning all the 
way from the disk to the event horizon). 

We should note that although the accretion rate over the course of our simulation is 
reduced by a factor of 4 contrasting the simulation with a/M — 0.998 to the one with 
a/M = 0, the evident lack of equilibrium in the former case prevents us from saying what 
will happen in the long run. It is possible that when the inner torus mass grows large enough, 
the time-averaged accretion rate will rise. In that case, the effect of black hole spin would 
be better described as altering the relation between accretion and inner disk mass rather 
than as a simple suppression of accretion. On the other hand, without doing more extended 
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simulations, we cannot say whether such an equihbrium is ever achieved — it is possible that 
the structure of the inner disk traverses some comphcated hmit cycle when the black hole 
spins very rapidly. In any case, it is clear that a magnetized spinning black hole can affect 
the rate at which matter accretes into it in ways that go beyond the changes in the effective 
potential usually considered. 

Another consequence of the substantial magnetic torques near the black hole is a re- 
duction in the net angular momentum brought to the hole as a result of accretion. As 
mentioned briefly in Paper 1 and at greater length in Gammie et al. (2004), the accreted 
angular momentum per unit rest-mass decreases so much with increasing a/M that the black 
hole can spin down (in the sense of decreasing a/M) even as M increases. Our results are in 
approximate agreement with those of Gammie et al. in indicating that there is near balance 
when a/M 0.9, but likely net spin-down for greater a/M. This result suggests that black 
holes cannot be spun-up by accretion to a/M > 0.9; only a black hole merger might create 
a more rapidly spinning black hole. 

We also examined energy transport in the vicinity of the black hole. As with angular 
momentum transport, the combination of magnetic fields and rotating holes leads to phe- 
nomena not normally included in standard disk theory. Magnetic fields can extract both 
energy and angular momentum from the black hole in a fashion resembling the Blandford 
& Znajek (1977) mechanism, but with several notable differences. First, the magnetic fleld 
is far from force-free. Second, the relevant fleld lines connect the plunging region to the 
disk and the outflow rather than only connecting the event horizon to inflnity (cf. similar 
results found by McKinney & Gammie 2004). Energy is drawn from a rotating black hole in 
a manner that, far from being independent of accretion as in the original Blandford-Znajek 
picture, is inherently dependent upon accretion. In fact, the most convenient way to describe 
its magnitude is relative to the rate at which rest-mass is accreted: for modest spins, the 
electromagnetic energy release per unit mass is < 1%, while for very rapid spin, it can be 
~ 10% or more. As the spin increases, the angular distribution of the Poynting flux also 
changes, with most going into the disk at low spin but the greater share entering the outflow 
when the black hole spins very rapidly. 

The simulations have demonstrated several new dynamical effects for magnetically 
driven accretion ffows in a Kerr metric. At this point we can give only qualitative assessment 
of the observational consequences of these effects. For example, strong outward magnetic 
angular momentum flux from a spinning black hole leads to a substantial alteration in the 
mass surface density distribution, creating a sizable inner torus in the region just outside 
rjns- Completely unexpected in conventional disk models, this large mass concentration 
could significantly alter the expected X-ray reflectivity of this region deep in the relativis- 
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tic potential. Because magnetic stresses continue throughout the inner disk and plunging 
region, the efficiency with which orbital shear energy can be dissipated into heat and then 
escaping photons could be greater than conventionally envisioned, but quantitative results 
require an improved treatment of thermodynamics. For rapidly spinning holes, the Poynting 
flux represents a very signiflcant part of the energy flow near the black hole. Although the 
impact of these effects on disk thermal emission remains to be worked out, their presence 
will surely lead to significant effects not found in traditional models. 
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A. Defining tiie Fluid Frame Tetrads 

In order to compute quantities measured by a physical observer, one must define a set 
of basis vectors (commonly called "tetrads") appropriate to that observer, which we write 
here as e^^^ for the contravariant set. In this notation, {u) identifies the direction in the 
observer's frame. 

Ideally, these should be directed in four-space so as to be most useful for physical 
interpretation. For example, in most instances it is important to be able to identify the 
direction of proper time in the observer's frame. This is simply done by setting e^^^ = C/^, 
the observer's four- velocity. In the observer's own frame, the only non-zero component of the 
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four- vector is the time component, so it automatically satisfies the requirement of appropriate 
directionality. Because U^U^ = —1, it is also automatically normalized to unit magnitude. 

Identifying other directions is, in general, more ambiguous. One can always construct 
the remaining three tetrads by a Gram-Schmidt procedure, requiring the newly-constructed 
tetrad element to be orthogonal to all the previously-defined ones and normalized to unit 
magnitude. However, there are many different ways in which the tetrad elements can be 
ordered for this procedure and the results depend on the ordering. In fact, the results depend 
in general on whether one first seeks contravariant or covariant versions. This mathematical 
ambiguity emphasizes the physical ambiguity in identifying spatial directions in one frame 
with those in another. 

For our purposes, we wish to construct a tetrad basis set that both reduces to the 
conventional directions in the Newtonian limit and may be most directly compared with 
the set chosen by Novikov & Thorne (1973). The choice made by Novikov and Thorne 
is a special case that is unambiguous because they were interested only in the limit that 
jjr — jjz _ Q the coordinate frame. When that is so, tetrad elements that are purely 
radial or purely vertical in the coordinate frame are automatically orthogonal to the four- 
velocity. Consequently, the R and z tetrad elements can be taken as purely radial and 
vertical in the coordinate frame, while the t and elements couple only to each other. This 
so simplifies the problem that no ambiguities remain. 

Although the components of the fluid four-velocities found in our simulations are in 
general non-zero everywhere, in much of the simulation region, U'^ ^ C/'', C/^. In an effort to 
ensure that the direction of greatest difference between the fluid frame and the coordinate 
frame is aligned as closely as possible, after setting e^^-^ = f/^, we begin the Gram-Schmidt 
process with the (f) clement of the tetrad and set its r- and ^^-components in the coordinate 
frame to zero everywhere. The r-clcment can then be constructed with a zero ^^-component, 
but it is also necessary to insist that e^^.^ > to preserve the right orientation of the tetrad 
set. Finally, only the ^-element remains to be computed; it has non-zero components in all 
four coordinates. The result is: 




[[/*, [/'■, U'f', U^] 



0, 



.g't"t>C2-g"^C,,0] 



(Al) 
(A2) 
(A3) 



-{1/N,) [Z, 0,1,0] 

(l/iVa) [l,iC^/Co){g''-l')/g"■,-l, 
C_/Cs - {U^/U'){C,/CoW - 1^)1 g'^] , 
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where 
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Co 
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Fig. 1. — The total mass found inside r = 15M as a function of time, given as fraction of 
the initial torus mass. The sohd hne corresponds to the a/M — 0.0 run, the dotted line to 
a/M = 0.5, the dashed hne to a/M = 0.9 and the dot-dashed hne to a/M = 0.998. In the 
three lower spin models the mass in the inner accretion disk remains roughly constant after 
t/M 2000, but the mass continues to grow throughout the simulation for the high spin 
model. 



-30- 



0.25 
0.20 

0.15 

0.1 

0.05 
0.00 







a/M = 0.0 



10 15 20 25 

r/M 



a/M = 0.9 




0.25 
0.20 
0.15 
0.10 

0.05 
0.00 







a/M = 0.5 



10 15 

r/M 



20 25 



a/M = 0.99f 




Fig. 2. — Time-averaged, shell- integrated accretion rate as a function of radius (solid curve), 
and rms accretion rate fluctuation about the mean at each radius (dashed curve). The units 
of accretion rate are fraction of the initial torus mass accreted per M of time (multiplied by 
10^). Plots for all four black hole spin models are shown. 
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Fig. 3. — (a) Time-averaged, density- weighted inflow velocity {{U^)p) as a function of r. The 
solid curve is the a/M = simulation, the dotted curve is a/M = 0.5, the dashed curve is 
a/M = 0.9, and the dash-dot curve is a/M = 0.998. (b) Ratio of inflow time (as defined 
by eqn. 6) to orbital period in the equatorial plane as a function of r/rms- The curves are 
identified as in (a). 
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Fig. 4. — Time-averaged surface density E as a function of proper radius, R, normalized to 
the marginally stable value, out to r = 15M in each of the four simulations. Each curve is 
also normalized to the initial torus mass of the simulation. The solid curve is a/M = 0, the 
dotted is a/M = 0.5, the dashed is a/M = 0.9, and the dot-dashed is a/M = 0.998. For 
comparison, the triple dot-dashed curve is the surface density predicted by a conventional 
Novikov-Thorne disk model for a/M = 0.9, assuming a constant ratio of stress to pressure. 
The amplitude of the curve is arbitrary; only the shape is relevant. 
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Fig. 5. — Absolute values of the shell-integrated and time-averaged angular momentum 
fluxes in the bound portion of the flow for the four simulations. The solid curve is the flux 
associated with the matter, the dashed curve is the magnetic torque, the dot-dash curve 
is the part due to advected magnetic energy. The sign convention chosen is such that the 
matter flux and advected magnetic energy carry prograde angular momentum inward, while 
the magnetic torque carries positive angular momentum outward. 
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Fig. 6. — The time-average flux of angular momentum per accreted rest-mass, {{T^))/{M) 
for each of the four runs. The data are time-averaged over the last 75% of the simulation. 
The dashed line is the angular momentum carried by the matter flux; the solid line is the 
total angular momentum flux. For both these curves the space integrals are restricted to 
the bound material, i.e., cells with —hUt < 1. The dot-dashed line shows the total angular 
momentum flux over the whole domain, including the jet. The values are contrasted with 
the angular momentum per rest-mass of the marginally stable orbit j(r^s) (dotted line). 
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Fig. 7. — The azimuthally-averaged stress as a function of polar angle 6, time-averaged 
over the last 75% of the simulation. The solid curves are the magnetic contribution to the 
outward-directed stress at r = 2rms, r = Vms and the inner boundary. In the a/M = 
model the magnetic stress is inward-directed (dot-dashed line). The dashed lines are the 
mass contribution to the stress, (i.e., phLf^U^), at r = r^s and the inner boundary. 
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Fig. 8. — The shell-integrated magnetic r-(J) stress in the fluid frame per unit area in the 
equatorial plane is plotted against radius for each of the four simulations at a selected late 
time. The dashed curves show the prediction of the Novikov-Thorne model for the time- 
average accretion rate at r^n for each of these simulations. 
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Fig. 9. — Specific energy fiux, defined as 1 — ((T/'))/(M) for each of the four simulations. 
The dashed fine is derived using the fluid energy flux only, while the solid line has both the 
fluid and electromagnetic energy flux. The dotted line corresponds to the speciflc energy 
of a circular orbit at each radius; the horizontal dotted line is the speciflc energy of the 
marginally stable circular orbit at at r^^. 



-37- 






5 10 15 20 25 



a/M = 0.9 




5 10 15 20 25 




0.000 0.002 0.004 0.006 0.008 0.010 



Fig. 10. — Time- and azimuthally-averaged outgoing Poynting flux in the four simulations, 
plotted on the same linear scale. The values are normalized by the initial energy and mul- 
tiplied by lO''. As the black hole spins more rapidly, the Poynting flux associated with the 
funnel-wall outflow rapidly grows in importance. 
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Fig. 11. — Ratio of time-averaged electromagnetic energy flux, {{T[))em to fluid energy 
flux, {(Tl)) fluid, as a function of radius. The solid line corresponds to the a/M = 0.0 run, 
the dotted line to a/M = 0.5, the dashed line to a/M = 0.9 and the dot-dashed line to 
a/M = 0.998. As the spin of the black hole increases, an increasing fraction of the energy 
flux is in the form of Poynting flux. 



